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ABSTRACT 

Aims. We study the structure of the magnetic elements in network-cell interiors. 

Methods. A quiet Sun area close to the disc centre was observed with the spectro-polarimeter of the Solar Optical Telescope on board 
the Hinode space mission, which yielded the best spatial resolution ever achieved in polarimetric data of the Fe 1 630 nm line pair For 
comparison and interpretation, we synthesize a similar data set from a three-dimensional magneto-hydrodynamic simulation. 
Results. We find several examples of magnetic elements, either roundish (tube) or elongated (sheet), which show a central area 
of negative Stokes-V area asymmetry framed or surrounded by a peripheral area with larger positive asymmetry. This pattern was 
predicted some eight years ago on the basis of numerical simulations. Here, we observationally confirm its existence for the first time. 
Conclusions. We gather convincing evidence that this pattern of Stokes- K area asymmetry is caused by the funnel-shaped boundary 
of magnetic elements that separates the flux concentration from the weak-field environment. On this basis, we conclude that electric 
current sheets induced by such magnetic boundary layers are common in the photosphere. 
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1. Introduction 

With ever increasing polarimetric sensitivity and spatial resolu- 
tion, it becomes now clear that even the most quiescent areas 
on the solar surface harbour ample amounts of magnetic flux. 
This flux becomes visible in small patches of field concentra- 
tions, called magnetic elements for short. While horizontal fields 
tend to occur near the edge of granules, the present investigation 
focusses on fields predominantly oriented in the vertical direc- 



tion and mainly occurring in the intergranular space (Lites et al. 
2007b a). There, magnetic elements are often visible in G-band 
filtergrams as point-like objects. In more active regions, mag- 
netic elements of extended size appear in the form of 'crinkles', 
'ribbon bands', 'flowers', etc. ( [Berger et al.|2004| l, which show 
a sub-structure consisting of a dark central area surrounded by a 
striation of bright material and further outside by a downdraft of 
plasma (Langangen et al. 2007 ). While Doppler measurements 
of this downdraft are at the limit of spatial resolution with the 
best ground-based solar telescopes, spectro-polarimetric mea- 
surements have less spatial resolution due to image degradation 
by atmospheric seeing over the required long integration times. 

Measurements with the spectro-polarimeter on board the 
Hinode satellite have a superior spatial resolution of approxi- 
mately 0.3" thanks to excellent pointing capabilities ( |Shimizu| 
|et al.||2()07] l and the absence of seeing. Using Hinode data, we 
show examples of unipolar magnetic elements of the network- 
ceU interiors that possess a distinct sub-structure, which is strik- 
ingly manifested in maps of the asymmetry of spectral lines in 
the circularly polarized light — the Stokes-V area asymmetry. 

The areas of the blue and the red lobes of Stokes-V profiles 
are identical when formed in a static atmosphere, but become 



asymmetric in the presence of gradients in the velocity and mag- 
netic field strength (jmln g et al. 1975 f [Solanki & Pa hlke 1988] 
San chez Almeida et aH1989( |Landolfi & Landi degl'Innocenti] 
1996|). Hence, the variation in the Stokes-V asymmetry across 
magnetic elements gives information on their magnetic field and 
plasma flow properties. Since this information is involved, we 
computed synthetic Stokes-V profiles of magnetic elements that 
occur in a three-dimensional simulation of magneto-convection 
and compare their asymmetry with the measured ones in order 
to understand their origin. 



2. Observation and data reduction 

The observations were caiTied out on 10 March, 2007 with the 
spectro-polaiimeter (SP) of the Solar Optical Telescope (SOT) 



on board of the Hinode space satellite (Kosugi et al. 2007 
[Tsuneta et al. 12007] Suematsu et al.|2007 l. The four Stokes pro 



files of the Fei630nm line pair were recorded in a quiet Sun 
area close to the disc centre. The integration time of 4.8 s led to 
an rms noise level in the polarization signal of cr ~ 1.2 x lO^^* I^. 
The data calibration for the SP is described in llchimoto et al.l 
(j2007j) and jTarbell ( |2006 l). The field of view comprises an area 
of 162" X 302" with a spatial sampling of 0.16" along the slit 
and 0.15" in the scanning direction. The spatial resolution of the 
resulting spectro-polarimetric maps is approximately 0.3", and 
the spectral sampling is 2.15 pm. The rms continuum contrast 
at 630.0 nm is around 7.5 %, far above the 3% value typically 
observed from the ground at the same wavelength (Rezaei et al. 
2007 ; KJiomenko et al. 2005 ). This same data set served Lites 
et al.| ( |2007a,b) in an investigation of the horizontal magnetic 
fields in the quiet Sun. 
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Fig. 1. Columns a-c: observational data obtained with the spectro-polarimeter of Hinode/SOT. Columns d andf: synthetic data from 
the 3-D MHD simulation. Columns e and g: same as d and f but after application of the SOT-PSF to the synthetic intensity maps. 
Rows from top to bottom: the sign of the Stokes-V area asymmetry for the line Fe 1 630.25 nm, where black and white correspond to 
negative and positive, respectively, the Fe 1 630. 15 nm line-wing velocity in km s"' , Vtot of Fe 1 630.25 nm in mA, and the continuum 
intensity at 630.0 nm. Contours refer to Vtot- Distance between tick marks is 0.5". 



3. Numerical simulation 

The three-dimensional magnetohydrodynamic simulation com- 
prises an area of 4 800 x 4 800 km^, corresponding to 6.6" x 6.6" 
and spans a height range of 2 800 km from the top layers of the 
convection zone to the middle chromosphere. Details of the sim- 
ulation, carried out with the CO^BOLD code, can be found in 



Schaffenberger et al. (2005 2006 1. 



Initially, a homogeneous vertical magnetic field with a 
strength of 1 mT was superposed on a model of relaxed ther- 
mal convection. After relaxation, fields of mixed polarity occur 
throughout the photosphere with an area-fractional imbalance 
typically of 3:1 for fields stronger than 1 mT. Similar polarity 
imbalances of even larger fields of view also occur in observa- 
tional data (Lites 2002) l. Because of the periodic side boundaries 
and the conditions that the magnetic field must remain vertical at 
the top and bottom boundaries, the horizontally averaged vertical 
net magnetic-flux density remains 1 mT throughout the simula- 
tion. This value compares very well with the mean flux density 
of 1.12 mT reported by Lites et aT] ( 2007a ) for the Hinode data 



set. The mean absolute horizontal field strength is 2.5 mT in the 
upper photosphere of the simulation, 2.5 times the net vertical 
flux density. |Lites et al. ( 2007a I report a ratio of 5. 



The grid constant of the computational domain in the hori- 
zontal direction is 0.055" (120^ grid cells), three times smaller 
than the spatial sampling of the Hinode SP data. The simu- 
lation was advanced for 4400 s of real time. The results pre- 
sented here refer to the last 2400 s, well after the initial relax- 



ation phase. Using the SIR code ( Ruiz Cobo & del Toro Iniesta 
1992 Bellot Rubio|2003) , we synthesized the Stokes profiles of 
the Fe 1 630 nm line pair from a large number of simulation snap- 
shots with a spectral sampling of 2 pm. Profiles were computed 
along vertical lines-of-sight for each grid point. To the result- 
ing intensity maps, we additionally applied a theoretical point- 



spread function (PSF) that models the spatial resolution capabil- 
ity of the SOT optics ( ,Wedemeyer-Bohm|2008| ). 

4. Data analysis 

The line parameters extracted from both the calibrated Stokes 
data of the observation and from the simulation include the area 
and amplitude asymmetry of Stokes-V , the line-core and line- 
wing velocities, and the total (signed) circular polarization 



(1) 



where Aq is the zero-crossing wavelength of a regular Stokes-V 
profile, and Ai- and A^ denote fixed wavelengths in the red and 
blue continua of the lines (similar to Rezaei et al. 2007). The 
relative Stokes-V area asymmetry ( |Solanki & StenflO|,1984| | is 
defined as 



6A 



= ( r ° I V(^)| dA- f \V(A)\ d/i) / r ' I V(^)| dA . (2) 

Figure [T] shows a few selected examples of magnetic ele- 
ments from the Hinode data (columns a-c) and the simulation 
(columns d-g). The observed ones are all taken from network 
cell interiors. Columns e and g are based on the same synthetic 
data as are columns d and f, respectively, but degraded with the 
theoretical PSF for SOT. The top row displays the sign of the 
area asymmetry: black corresponds to negative and white to pos- 
itive area asymmetries. Gray indicates that the Stokes-V profile 
has a signal below the 3cr noise level (65 % of the total field of 
view) or that it has an irregular shape (20 %), where most of the 
latter would classify as regular when lowering the noise level. 
This same threshold also applies to the simulation data of the 
panels in columns d-g, where we ignore V-profiles with an am- 
plitude less than 3.6 x 10"^ /c for comparison with the Hinode 
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Fig. 2. Variation in 6A across magnetic elements from the 
Hinode data (top row) and the simulation (bottom row). From 
left to right and top to bottom, the sections are taken from the 
data corresponding to columns a, b, d, and f of Fig.[T| 



data. In the rightmost two panels, this threshold is 5 times lower 
in order to highlight the case of a weak magnetic element that 
would have barely been detected with Hinode. 

The line-wing velocity of Fei 630.15 nm, shown in the sec- 
ond row, is obtained by averaging the Stokes-/ Doppler shift 
in the intensity interval 0.8-0.9/c. We set the convective blue 
shift of the line-core velocity in the quiet Sun to -0.2 km sec 
Because of the low degree of polarization, the line-wing velocity 
is not spoilt by magnetic influence. 

The third row displays Vtot, whose contours are shown in all 
the panels. It is a linear measure of the magnetic flux density 
according to Lites et al. ( 1999[ l and Lites et al. ( 2007b[ l. From 
the fourth row, which shows the continuum intensity at 630 nm, 
we see that the magnetic flux concentration occurs mainly in the 
intergranular space. 



5. Results and discussion 

Both the observation and the simulation show magnetic elements 
of either sheet like or roundish shape, possessing the same strik- 
ing pattern with regard to the Stokes-V asymmetry: a central re- 
gion with negative Stokes-V area asymmetry is framed or sur- 
rounded by a narrow seam of pixels having positive area asym- 
metry (Fig.[T] first row). Plotting 6A across a magnetic element, 
we find that |&4| is typically larger in the periphery than in the 
central region of the magnetic flux concentration. Representative 
cross sections that demonstrate this behaviour are shown in 
Fig. |2] for both the observation and the simulation. Along with 
6A in the top right panel of Fig. |2] Fig. 4] also gives the corre- 
sponding observed Stokes-V profiles. A so from Fig. [T]we see 
that the magnetic elements are located within and suiTounded by 
a downdraft. 

For gaining deeper insight, we now consider in Fig. [3] the 
vertical cross section through the simulation snapshot shown in 
column d of Fig. [T]for the position y - 1.15". It displays the 
magnetic field strength (colour-coded) together with the veloc- 
ity field, projected on the vertical plane (white arrows). The hor- 
izontally running curve indicates optical depth r^oonm = 1. A 
magnetic flux concentration is forming in the downdraft and has 
reached a field strength of 0.13 T at rsoonm = 1- Close to the sur- 
face of optical depth unity and through the photosphere, a rela- 
tively sharp, funnel-shaped boundary separates the magnetic flux 
concentration from the weak-field or field-free surroundings. 
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Fig. 3. Vertical cross section through the simulation box cor- 
responding to position y - 1.15" of the column d of Fig. [T] 
It displays the logarithmic magnetic field strength (from 1.0 to 
10^^ mT, colour-coded), together with the velocity field, pro- 
jected on the vertical plane (white aiTows) and the electric cur- 
rent density normal to the plane (black contours). Optical depth 
T"500nm = 1 is shown as well. The white vertical lines indicate 
the ranges in which Stokes-V profiles from vertical lines-of-sight 
have either positive or negative area asymmetry, 6A. Outside this 
range, the Stokes-V signal is below the 3cr noise level of the ob- 
servations. dA{x) is plotted in the lower left panel of Fig. |2] 



The inner vertical white lines indicate a central part of the 
flux concentration (880 < x < 1120 km), where Stokes-V pro- 
files have 6A < 0. This is because their Unes-of-sight sam- 
ple increasing field strength and increasing downflow veloc- 
ity with increasing continuum optical depth ( So lanki & PahUEel 
1988 ). This situation changes drastically in the peripheral region 
840 < X < 880 km and 1120 < jc < 1280 km (also indicated 
with white lines), where M > 0. Here, the downflow speed still 
increases with optical depth; but because of the funnel-shaped 
flux concentration, lines-of-sight traverse a magnetic boundary 
layer, where the field strength rapidly drops with increasing op- 
tical depth. The combination of increasing velocity and decreas- 
ing field strength produces positive area asymmetry. From this, 
we conclude that the narrow seam of positive area asymmetry 
in the periphery of magnetic flux elements seen in the Hinode 
data is a consequence of their well-defined boundary (or magne- 
topause) that expands with height in a funnel-like manner. 

As an immediate consequence of Ampere's law, this bound- 
ary carries an electric current sheet. The black contours in Fig. |3] 
indicate a current density of 0.5 Am"^, encircling higher values 
of up to 4.0 Am"^. This refers to the current component perpen- 
dicular to the section plane, which flows in the opposite direction 
on either side of the magnetic flux concentration. 

If the magnetic flux concentration is inclined enough with re- 
spect to the viewing direction, lines-of-sight traverse the bound- 
ary layer on only one side of the flux concentration, which leads 
to a one-sided seam of positive area asymmetry. This is almost 
the case in Fig.[3](and corresponding M-panel of Fig. [T}, where 
the flux concentration is slightly inclined to the right making the 
seam of positive 6A to the left much narrower than that on the 



Rezaei et al.: Variation of area asymmetry across a magnetic element 





0.10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.26 O.JO 0.10 0.15 0.20 0.25 0.30 0.10 0.16 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 
[\ - 6301 (nm) [^ - 6301 (fim) \X - 5301 (ntn) f\ - 6301 [nm) f^ - 6301 (nm) fX - 530] (nm) f^ - 6301 (nm) 

Fig. 4. Stokes-V profiles from the Hinode data across the magnetic element of column b of Fig. [l] 6A(x) for the same section is 
shown in Fig.|2](top right panel). 



right side. In fact, most magnetic elements in the field of view of 
the observation show a one-sided pattern like in the upper right 
corner of the synthetic M-map in columns d or e. 

Because the mean magnetic flux density in the network cell 
interior is low, granular flow gathers only small amounts of 
flux so that only small and weak magnetic flux concentrations 
form dSteiner 2003 1. They form and disperse with the onset 



and cessation of intergranular downdrafts. Therefore we observe 
strong downflows associated with the magnetic elements con- 
sidered in this work. Indeed, Grossm ann-Doert h et al.| ( |1996| l, 
Sigwarth et"ar ('1999'), and Martinez Piflet et al. ( 1997 ) also ob- 
serve mainly downflows associated with weak magnetic signals. 
In combination with increasing field strength, these downflows 
give rise to negative area asymmetry. Yet, the peripheral posi- 
tive area asymmetry is larger (Fig.[2]i and may partially balance 
or even outweigh the negative contribution when observing at a 
lower spatial resolution. 

Langangen et al. ( 2007| l finds that the downdraft velocity 



near the edges of magnetic elements is larger than in the cen- 
tral part. Here we find in both observations and simulation the 
maximum speed in the centre of the flux concentration. In the 
simulation this is because the magnetic element forms as a con- 
sequence of an intergranular downdraft, with peak velocity in 
the centre. This is diff'erent from the regime of a mature strong 
flux concentration for which simulations exhibit veritable down- 
flow jets near their edges but not in the centre ( [Steiner et al.| 
1998 Shelyag et al. 2007 1. This regime might be more represen- 



tative of the active region elements observed by [Langangen et aL 
( |2007ll. 



" Grossmann-Doerth et al.| ( fT988l [1989) 1 and |Solanki| ( fT989l l 
first pointed to the peripheral, canopy-hke magnetopause of 
magnetic elements as the origin of the observed positive Stokes- 
V area asymmetry. [Sterner (fT999| found this effect at work in a 
two-dimensional MHD-simulation, but also found that the cen- 
tral region of the magnetic flux sheet tends to have negative val- 
ues and that a delicate balance with the peripheral region ex- 
ists in which the positive values may dominate and outweigh 
the negative ones. These findings have recently been confirmed 
in a three-dimensional environment by S helyag et al.| (2007) . 
[Bellot Rubio et al.| ( |2000j came to a similar conclusion based on 
semi-empirical modelling. [Leka & Steiner ( |2001[ ) found a simi- 
lar pattern in the Stokes-V area asymmetry in pores and magnetic 
knots. Here, for the first time, we find this pattern in observations 
of small-scale magnetic elements in the network cell interior. 



6. Conclusion 

Spectro-polarimetric data of a quiet Sun area close to the disc 
centre obtained with SOT on board Hinode were analysed. We 
find magnetic elements of either a roundish (tube) or an elon- 
gated (sheet) shape, which show Stokes-V profiles of negative 
area asymmetry in the centre, surrounded or framed by pixels 
of larger, positive area asymmetry. A comparison with results 



from 3-D MHD-simulations suggests that this peculiar pattern 
in Stokes-V area asymmetry (first predicted by Steiner 1999) is 
due to the confined nature of the field in magnetic elements with 
a funnel-shaped boundary layer that gives rise to a steep gradient 
in field strength along the line of sight. 

We also conclude that this kind of magnetic element of the 
internetwork is accompanied by electric current sheets. While 
these conclusions are evident from the comparison, we cannot 
exclude the hypothesis that a suitable magnetic structuring on 
scales not resolved by the present observation and simulation 
would lead to the observed pattern in 6 A. 
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